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Abstract 

The sinh-Gordon model is restricted to a half-line by boundary conditions main- 
taining integrability. A perturbative calculation of the reflection factor is given 
to one loop order in the bulk coupling and to first order in the difference of the 
two parameters introduced at the boundary, providing a further verification of 
Ghoshal's formula. The calculation is consistent with a conjecture for the gen- 
eral dependence of the reflection factor on the boundary parameters and the bulk 
coupling. 
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1 Introduction 



Over the last few years, following the pioneering ideas of Ghoshal and Zamolodchikov, 
and others |l|-|5[], much work has been done to investigate integrable quantum field 
theory with a boundary. In particular, the affine Toda field theories have offered a 
surprisingly rich structure which is just beginning to be understood. The classical affine 
Toda field theories are known to remain integrable in the presence of certain (generally 
quite restricted) boundary conditions confining them to a half-line, or to an interval 
P-TT5I]. The corresponding quantum field theories are hardly explored although there 
has been progress in certain cases associated with the an class of models P, |l6|-|i"8f . 
The first of these (n = 1) is the sinh-Gordon model and, unlike all the others has a 
set of integrable boundary conditions depending on a pair of extra parameters, the so- 
called boundary parameters. However, even in this case, it remains to be seen precisely 
how the two boundary parameters influence quantities of interest such as the reflection 
factor. 

The sinh-Gordon model is related to the sine-Gordon model and its reflection fac- 
tors are related to the reflection factors of the lightest breather in the sine-Gordon 
model. Indeed, the results obtained by Ghoshal and Zamolodchikov [l|,|2|] allow a de- 
termination of the general form of the sinh-Gordon reflection factor and this is a very 
useful piece of information. However, apart from two special cases (Neumann and 
Dirichlet boundary conditions) Ghoshal and Zamolodchikov's formulae fail to provide 
a relationship between the reflection factors and the boundary parameters themselves 
which appear as the data in a Lagrangian formulation of the model. 

Recently, it was noticed [19| that for certain ranges of the boundary parameters in 
the sinh-Gordon model there are real periodic classical finite-energy solutions called 
boundary breathers. The sinh-Gordon model has no real finite-energy solutions at all 
on the whole line (other than = 0) but, once there is a boundary singularities may 
be hidden behind it allowing solutions restricted to a half-line to have finite energy. 
The existence of periodic solutions allows a semi-classical quantization determining 
the spectrum of boundary bound states. Once the spectrum of bound state energies 
is known it may be compared with a boundary bootstrap calculation of the energies 
of the same states and hence may be used, in principle, to find a relationship between 
the Lagrangian boundary parameters and the data in Ghoshal's formula. In fact this 
is complicated and was carried out in [19] for the simpler situation in which the two 
boundary parameters are equal and the —>■ —0 symmetry of the bulk theory is pre- 
served. The corresponding calculation in the more general case is not yet completed. 
The semi-classical calculation may turn out to be exact (as it did for the energy spec- 
trum of the breathers in the bulk sine-Gordon model itself; for a review see |2(J) but, 
failing a proof of exactness the results may be checked against low-order perturbation 
theory 

If the two boundary parameters are different and the — > — symmetry is broken, 
the perturbation theory becomes substantially more complicated: the lowest energy 
static background is no longer the configuration = and therefore the perturbation 
theory must be developed within a non-trivial background, leading to additional cubic 
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and higher odd order couplings as well as a substantially more intricate propagator. 
Nevertheless, as will be shown below, provided calculations are restricted to first order 
in the difference of the two boundary parameters, some of the complications disappear 
and we are able to calculate the correction to the reflection factor at one loop. The 
result we obtain is consistent with a conjecture for the relationship between the bound- 
ary parameters and Ghoshal's formula in the more general setting. A full discussion, 
to all orders in the difference of the boundary parameters, even at one loop order in 
the bulk coupling, must be postponed. 



2 sinh-Gordon model 

The sinh-Gordon theory corresponds to the affine Toda field theory associated with 
The bulk Lagrangian density of the theory is: 

C^\d^<j>-V{<j>) (2.1) 

where 

2m 2 

V(<P) = ^cosHp a( p) (2.2) 

The real constants m and f3 provide a mass scale and a coupling constant, respectively, 
and it is customary in affine Toda field theory to choose a = 

Assuming two sinh-Gordon particles of relative rapidity G scatter from each other 
elastically, the bulk S-matrix characterizing this process is given by [p5| -27 



s < e > = -jByhw) < 2 ' 3 > 

where 

1 3 2 

_ i_ (2 4.) 

27rl + /5 2 /4vr 1 ' ' 

and the symbol ( ) denotes the hyperbolic building block: 

_ sinh(6/2 + g) 
1 j "sinh(e/2-^f) 1 j 

The S-matrix is invariant under the following weak-strong coupling transformation 

3 -> Att/(3 (2.6) 

a property known as weak-strong coupling duality. 

A sinh-Gordon theory on the half-line is described by the following Lagrangian 
density: 

C = 9(-x)C-8(x)B, (2.7) 
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where B is regarded as a functional of the field only, not its time derivative. Moreover, 
the generic form of B is given by Jl| 

6 = ^L e -^ + (7ie ^) (2.8) 



ft 2 ■ 

Note, the coefficients <To and o"i are a pair of real numbers, essentially free (but see 
J7, 1C]), which represent the extra parameters permitted at the boundary x = 0. 

The equation of motion and the boundary condition for the sinh-Gordon model on 
the half-line become (after rescaling the mass scale to unity): 

^ = _^^*_ e -^ x<0 (2.9) 

^ = -^( aie WV2_ aoe -WV2\ X = Q (210) 

ox ft V / 

3 Reflection factor 

When a sinh-Gordon particle is moving towards the boundary located at x — it will 
reflect elastically from it meaning that the in- and out- one particle states, conveniently 
labelled by rapidity, will be related by a reflection factor 

\6> out =K(6)\-6> iQ . (3.1) 

The general form of this reflection factor is known although its detailed dependence 
on the boundary parameters <To and o\ appearing in (|2.8|) is known only in special cases. 
The sequence of arguments determining its form is somewhat indirect, stemming from 
the work of Ghoshal and Zamolodchikov as follows. Solving the boundary Yang- 
Baxter equation, and using general constraints implementing unitarity and a form of 
crossing symmetry, it proved possible to calculate the reflection factor for the sine- 
Gordon soliton. Since breathers are soliton-anti-soliton bound states, a subsequent set 
of calculations using the boundary bootstrap led Ghoshal [[J to conjecture reflection 
factors for each member of the full tower of breathers. Finally, the reflection factor for 
the lightest breather is supposed to be identical with that of the sinh-Gordon particle 
provided the sine-Gordon coupling ft is replaced by ift. 

Thus, suitably transformed in the manner described, Ghoshal's formula for the 
sinh-Gordon reflection factor is given by: 

K (e) = (l)(2-E/2)(l + g/2) 

g{ ] (l-%(7 1) i3))(l + %ff^))(l-FK(r 1 ^)(l + FK(7 1 ^) 1 ' ' 

Actually, Ghoshal's notation was a little different and made use of two other quantities 
rj, *d defined by E = Br]/n,F = iBd/it. 

There are special cases in which E and F have been conjectured. For example, the 
Neumann boundary condition which is defined by 

^ = when x = 0, (3.3) 
ox 



has been argued by Ghoshal and Zamolodchikov to need a reflection factor containing 

F = 0, E=l-B/2. (3.4) 



More recently |T9 |, the boundary breather states of the sinh-Gordon model re- 
stricted to a half-line were investigated and their energy spectrum calculated in two 
ways. First, by using the bootstrap equations, and then by finding a set of periodic 
finite-energy solutions which could be quantized using a WKB approximation. In- 
sisting that the two methods agreed with each other provided strong evidence for a 
relationship between the boundary parameters, the bulk coupling constant, and the 
parameters appearing in the reflection factor. For technical reasons this work was re- 
stricted to the case 00 = (J\ but, nevertheless, yielded an expression for E, with F = 0. 
Specifically, setting a = <j\ = cos an and restricting a to the range 1 > a > 1/2, E is 
given by, 

E = 2a(l-B/2). (3.5) 

In the limit a — ► 1/2 from above, this is in agreement with ( |3.4T ). 

Both of these conjectures are underpinned by low order perturbation theory. Kim 
2~Tf has calculated the one loop correction to the classical Neumann boundary condition 



reflection factor and found agreement with Ghoshal's formula to 0(f3 2 ). On the other 
hand, a more general calculation was carried out, also to one loop, agreeing with the 



conjecture (|3.5|), and indeed, preceding it |23]| . In fact, the perturbative calculation 
in [23] agrees with ( |3.5|) for a in the range 1 > a > 0. 

However, when a ^ ai the situation becomes much more complicated. Indeed, all 
that is known up to now is the behaviour of E and F in the limit /3 — > 0, deduced from 
a direct calculation of the classical reflection factor in the general case J7|. To describe 
this limit it is convenient to set cr = cosa 7r and <j\ = cosa!7r, with |a$| < 1, i — 0, 1. 
Then, 

£ , (0,cr ,(Ti) = a + ai, F(0, a , cxi) = a - ax. (3.6) 



and, the classical reflection factor itself is given in term of the basic factors (|2.5|) by 

(l) 2 

KJB) = —. hr, 77 7- ( 3 - 7 ) 

(1 — a — a x )(l + a + ai)(l — a + ai)(l + a — ax) 



The expression ( p.7|) is an essential ingredient of the basic two-point function, or 
propagator, which takes the form 



, . dudk i e -Mt-t') 




x (f(k, x)e ik{x - x,) + K c f(-k, x)e- ik[x+x,) ) (3. 



where 



ik + 2 K J 



The parameter x enters the static background 



e M V2 = £ + ^ (31Q) 

1 _ p 2(x-x ) v > 
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and determines the point at which the background becomes singular. For this reason, 
it is crucial that x > 0. Actually, x is determined by the boundary condition ( [2.10| ) 
and satisfies 



cothxo 



l + (7i 



(3.11) 



On the understanding that o\ > <jq, it is enough to consider this situation since ( 2.1U| ) 
is symmetric under 0^—0 and the interchange of cq with o\. 

In the expression ( pTfl ) for the propagator, the classical reflection factor appears as 
the coefficient of the reflection part of the free field two-point function calculated within 
the classical static background. To check Ghoshal's formula, the strategy introduced by 
Kim |21[ and developed in |23[ is to calculate perturbative corrections to the two-point 
function and then to identify corrections to the classical reflection factor by picking 
out the coefficient of e~ lk ( x+x ) as x,x' — > — oo. 



4 First order quantum corrections to the reflection 
factor 

Previously, as noted above, the first order perturbation calculation was achieved for 
the specially symmetric case (Tq = a%. In this section, this restriction will be relaxed 
and enough of the ingredients of the perturbation expansion will be calculated to allow 
a calculation (in principle) up to (3 2 . We will discover that generally the theory needs 
three- and four- point couplings which depend upon the x-dependent static background. 

Expanding the bulk potential ( |2T2| ) around the background solution to the equation 
of motion, (po(x), we derive the three- and four-point couplings: 



and 



^i = ^smh(v^0o), (4.1) 

^l = ^ 2 cosh(v / 2/30o)- (4.2) 

On the other hand, the static background is ( |3.10| ). So, after some manipulation, 
we have 

Cfik = cosh 2{x - x ) (coth 2 2{x - x ) - l) , (4.3) 

and 

C£l k = ^ 2 (2coth 2 2(x-x )-l). (4.4) 



In the same manner, using (|2.8| ) we can derive the three- and four-point couplings 
associated with the boundary term: 

,(3) _ V%P_ ( pl 3^ /V2 _ -„-!3<h/V2\ 



= ^ {^ 0/V2 ~ W-W) , (4.5) 

t f aie ^0/V2 + aQe -^o/V2\ (4 6) 

48 V / 



and 

boundary 
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For future reference, it will be useful to know these to first order in the difference 
of the two boundary parameters e = er — <j\. 



and similarly, 



cffi, = ^7+"^' + - (4-7) 

3 1 + 01 

^1 = ^ 2 + - (4-8) 
c ® _ ( <_) , u o) 

^boundary ~ u ^ 1+(J J +- ^ J > 
^^ = ^(2^1 +6) + .- (4.10) 

We shall also need the expressions for f(k, x) (|3.9| ) and the classical reflection factor 
K c ( |3.7| ) to the same order in e since both contribute to the propagator (|3.8|). In fact, 

f(k,x) = l + 0(e 2 ) (4.11) 

and the classical reflection factor ( |3.7| ) reduces to 

ik + 20 2ik ^ , 9 . . 

K c =- + -^ r^e + Oe 2 . 4.12 

ik-2a (ik-2a) 2 V ; V ; 

It is also convenient to write 

K c = K + eK h (4.13) 

where K is the classical reflection factor when the two boundary parameters are equal. 

To calculate quantum corrections to the classical reflection factor at one loop order 
we use perturbative methods generalised to the affme Toda field theory on a half-line 
21|-p4|. (For earlier references on boundary perturbation theory in general see; |28|-|30] 



for affine Toda perturbation theory see , or the review HO The 0{f3 2 ) correction 



to K has been calculated before and the purpose of this article is to calculate the 
corrections to K\ to the same order. 
The possible diagrams to 0(j3 2 ) are: 



I II III 

Figure 1: Three basic Feynman diagrams in one loop order. 
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These will be computed in configuration space noting that each vertex may either be 
situated at the boundary or within the bulk. In effect, there are ten contributions 
which need to be calculated. 

However, there is a simplifying feature provided we are content to work to first 
order in e. To recognise this it is enough to note that because the three-point coupling 
is already 0(e), implying the type II and III diagrams involve e 2 , only diagrams of type 
I need concern us. 

Thus, there are two contributions which need to be calculated in order to be able 
to deduce the quantum corrections to the classical reflection factor, both of type I. The 
first is directly related to the boundary, when the interaction vertex coincides with the 
boundary, and it takes the form 

ao2 r+oo 

-^(2<n + e) / dt"G(x,t;0,t")G(0,t";0,t"),G(0,t";x',t'), (4.14) 

including the correct coupling constant and combinatorial factors. 

The second contribution refers to the bulk potential which means the interaction 
vertex is located inside the bulk region x < 0. This contribution is given by 

p+oo pO 

— 4i(3 2 / dt" / dx"G(x, t; x" , t")G{x" , t"; x" , t")G(x" , t"; x' , t'). (4.15) 



Again, the combinatorial factor has been included. 

Let us first calculate the boundary contribution ( [4.14| ). The loop propagator is 
given to 0(e) by 



Gmn ^JJ^^X^ + ^ + ¥ ^), (4,6) 
where we have defined 



P (u", k") = — — . (4.17) 

uv ' ; u" 2 - k" 2 -4 + ip v ; 

Note, the integral is clearly divergent but the divergence is removed by a renor- 
malization of the boundary parameter. In effect, rearranging the part of the integrand 
containing the offending terms as follows, 

ik"-2a x ik" -2o x v ; 

and making a minimal subtraction of the divergent portion, (i.e. deleting the '2'), 
renders the integral finite. 

The part of the integral independent of e has been calculated before , therefore 
we may write, 

„, „ „. aicosai7r f f duo" dk!' , „ , „, 2ik"e 

G(0 - t iat > = -kzr + 'J J ^ p ^- k) w^- (419) 
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and concentrate on the 0(e) part. 

As far as the remaining integral is concerned, focusing on the energy variable and 
closing the integration contour in the upper half-plane, we encounter a simple pole at 
V k" 2 + 4. Therefore, we need to evaluate 

1 fdk" 1 2ik" 

2 J ^TTPri (ik« -2o x y € - ( - 20) 

The k" integration may be performed by closing the contour in the upper (lower) half- 
plane depending on whether o\ > (eri < 0). However, because of the branch points 
at ±2z it is convenient to locate the associated branch cuts along the imaginary axis 
from ±2z to oo. Evaluating the integral and assembling the pieces, we obtain 

duj" dk" „ . „ , 2ik" e ai e cos ai7r 1 , , , 

P (^ , k")— — e = ^ — + = . (4.21) 

2tt 2tt uv ; (iA; // -2a 1 ) 2 2sin 3 ai7r 2vr sin 2 ai7r V ; 

At this stage, the contribution to ( (4.14 ) is 

f3 2 ( aicosai7r e a\ ecosai7r 
-»— (2<ri + e) ( — — — - - + — —— 2 




4 \ sinai7r 2 sin 2tt sin ai7r 



lib - 2ax) 2 

The integration over i" ensures energy conservation at the interaction vertex and creates 
a Dirac delta function which immediately removes one of the energy variables, for 
example u' . The remaining integral over the momenta k and k! can be performed by 
completing the contours in the upper half-plane and taking into account the poles at 
k = k' = \J iJ 1 — 4 = k. However, if a\ > it is evident that the expressions for 
Kq and K\ have no poles inside the contour. If <j\ < 0, there is an additional pole 
but its contribution turns out to be exponentially decreasing in the asymptotic region 
x, x' — ► — oo. 



Finally, we obtain the boundary contribution Q4.14 ) in the form 

_ { f?_ f duj c _ iuj(t _ t , )c _ ik(x+xl) { 2a 1 cos 2 a 1 7r 1 

4 J 2tt \ sinam {ik-2^) 2 

aiCOsaiiT aicosaiir cos 2 ai7r \ e 



sinai7r sin 3 ai7r irsm 2 aiT[J (ik — 2ai) 2 



+ 4aiC ° s2Ql7r , £ }, (4.23) 

sinaivr (jfc _ 2ai) ? 



where k = 2 sinh 0. 

Next, we need to calculate the contribution (|4.15|) which to 0(e) is: 
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-4t(3 2 I dt" I dx" I I ^^ e -***-^iP (u>,k) 



CO 






2tt 2tt 

e -ik(x-x") | ^ + c -ik(x+x") , 2z£;e e -ifc( a; + a ;") 
iJfe - 2(Ti (iife - 2m) 2 

^^p (o/', k") ( i + ^ + 2 %-^y + aFc ^ 

2tt 2tt ov ' ; V iife" -2m (ifc"-2m) 2 



2tt 2tt 



iP (u', k') e~ w 



6 + tk'-2a 1 6 + (tk'-2a i y 6 J { ' 

The integral over t" yields a delta function which replaces u' by u. Furthermore, to 
calculate the integration over x", it is convenient to use the following device 

° dx"e ikx " +TX " = (4.25) 

where the small positive quantity r will be taken to zero at the final stage of the 
calculations. 

The loop integral which corresponds to the middle propagator of (|4.15| ), is obviously 



logarithmically divergent and this divergence will be removed by the infinite renormal- 
ization of the mass parameter in the bulk potential. So, after making the subtraction 
and integrating x" and u", and as before concentrating on the 0(e) piece, we obtain 
the contribution 



/ 



i_rr rdudkdk!_ {t _ tl) i{kx+k , xl) k)lPo{u k >) 

2 J J J 2tt 2tt 2tt V ' ; V ' ; 

dk" 1 f 2ike ik" + 2a x ( 1 



+ 



+ 



2tt ^k" 2 + 4 I (ik - 2o x ) 2 ik" - 2m \-k + k' - 2k" - ir 

1 ik' + 2cn 

~k + k' + 2k" + ir ik' - 2o x 

2ik'e ik" + 2a x ( 1 1 ik + 2a x 



[ik' - 2m) 2 ik"-2a x \-k' + k - 2k" - ir k' + k + 2k" + ir ik-2a x 
2ik"e ( 1 1 ik! + 2a x 



[ik" - 2m) 2 V k + k! - 2k" - ir k - k' - 2k" - it ik' - 2a x 

1 ik + 2a x 1 ik + 2a x ik' + 2m 



k-k' + 2k" + ir ik - 2o x k + k' + 2k" + ir ik - 2m ik! - 2a x 
In order to evaluate the integral over k", we encounter the following two types: 
dk" 1 / ik" + 2a x \ 1 



(4.26) 



2tt ^k" 2 + 4 \ik"-2a x ) (k + k' - 2k" - ir 



(4.27) 



and 

dk" 1 2ik"e 1 

(4.28 



2tt ^/k" 2 + 4 (ik"-2a 1 ) 2 (k + k' - 2k" - ir) 

Both of these may be performed by closing an appropriate contour in the upper half- 
plane ensuring that it runs around the branch cut located from k" = 2i to infinity 
along the imaginary axis. Note, if o\ > 0, then there is no pole inside the contour; 
however, if a\ > 0, there is an extra pole but its residue integrated over k and k! will 
give a contribution vanishing in the limit x, x' — > — oo. 
For example, the integral ( |4.28|) evaluates to 

dk" 1 2ik"e 1 e k + k' a l7 i 



2tt ^/k" 2 + 4 (ik" - 2a 1 ) 2 (k + k' - 2k" - ir) ir (ik + ik' - 4cxi) 2 sin 
e io\ 1 e ia\ aiir 

7r (40"! — ik — ik') sin 2 a\ix n (Aoi — ik — ik') sin 3 a X TX 



v t i ln Wgjjg^ , (4 , 9) 

7r (zA; + iA;' - 4(T!) 2 _ / (k+k') 2 | ^ ) 1 + ' (fc+fc,) - i / ( fc + fc 2 | 4 



Let us divide the bulk contribution ( f4.26| ) in two parts: -Bi containing integrals 



over /c" of type (|4.27| ) and B2 whose k" integrations are of type (|4.28|) . For both it is 



necessary, after performing the integration over cu", to do the k and k! integrals via 
closing contours in the upper half-plane to pick up the poles at k or k' — ^/u 2 — 4 = k. 
All other pole contributions lead to exponentially damped terms in the limit x, x' — > 
—00. 

After some manipulation, Bi is found to be equal to 

B x = -2f3 2 [ — c ~i"(t-t') c -ik(x+x') 1 2lkt 



2 



2vr (2A;) 2 (iifc - 2^; 

* ' mi %k +~ ^ (--^[. (4.30) 




4 sinai7T ifc - 20"! vr y/p + 4 V 2 



Notice that in the last term inside the braces, depends on 9 in a manner which 
potentially is very inconvenient for a comparison with Ghoshal's formula. Fortunately, 
this term will be canceled by a matching term in B2. 

After somewhat lengthier calculations, B2 is given by 

du ■ ,,, , „/\ 1 —2ie a 2 



B 2 = -2(3 2 [ ^ e -Mt-t') e -ik(x+*')_ 

J 2n 1 '- » « -™ 2 



(2A;) 2 I (iifc - 2(Ti) 2 7rsm z ai7r 



2ze diof 2iefc 



k - 2a x ) 2 sin 3 aiTr (ik - 2ai) 2 tta/p + 4 

+ i±^-( % 1 + ieaiai ))■ (4.31) 

ik — 2o\ \ 27rsin 2 ai7r 2 sin ai7r > 
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Assembling the 0(e) part of flOg ) with ( P~30l) and ( P~3~T1 ) we obtain, 



+ 



2 ./ 2tt 
1 



g -iuj(t—t')^—ik(x+x') 



d\0\ 



k — 2e>"i) 2 \ 27r 2sinai 



7T 



(— 2ai sin ai7r) 



, (4-32) 



(z/c - 2a!) 3 x ' ' ■ (ife-2 < 7 1 ) 2 fe\ 4 2^ + 4, 

whence we can deduce the correction to the quantity K 1 in ( |4.13| ). Explicitly, we have, 

1 ( 1 Oiffi \ 1 



5K t 



-i(3 z ek 



I (XiOi 

+ „ . + 



{ik-2a l ) 2 \^ 2sina 1 7ry (ik-2a x f 
1 'lit i 



(— 2ai sin ai7r) 



+ 



2 / i 



ik-2a 1 ) 2 k \ 4 + 4 



The correction to Kq which was calculated before in |23 

1 



6K 



id 2 

-—K {k) sinh# 



+2tt! 



is. 
1 



cosh 9 + 1 cosh 9 
1 I 



cosh 9 — sinaiTT cosh 9 + sin 
This completes the collection of ingredients we need. 



(4.33) 



(4.34) 



5 Comparison with Ghoshal's formula 

In this section, the corrections to the classical reflection factor calculated above will be 
compared with the formula of Ghoshal quoted in Q3.2|) . 

Using Q4.13 ), the relative correction to the classical reflection factor K c is given in 
terms of the corrections SKq and 5K\ by 



SK C 



K l 5K Q + e (K^S^ - KiKq 2 5K ) . 



Hence, using Q4.33|) and (|4.34j) we have 



SKc 



- sinh 9 



-2a x 



cosh 9 + 1 cosh 9 
1 



1 



cosh 9 — sin a\n cosh 9 + sin a\i\ 
i/? 2 esinh# (1/1 1 



lsinai7r [ 7T \ cosh 6* — sina^ cosh 6* + sina^ 

1 1 

a\ cosai7r 



(cosh 9 — sin aiix) 2 (cosh# + sinai7r) 2 



(5.1) 



(5.2) 
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On the other hand, Ghoshal's formula ( |3.2| ) for the reflection factor up to one loop 
order is given by: 



K q {9) ~ K c {9) (l - ^ sinh 9 F(9)^ , (5.3) 



where 

HO) 



1 i 



cosh 9 + 1 cosh 9 

ei ei 



cosh^ + sin(e 7r/2) cosh# — sin(e 7r/2) 

I fi f i (54) 

cosh# + sin(/o7r/2) cosh^ — sin(/ 7r/2) 

In calculating ( |5.4j) we have made use of the expansions of -E 1 and F to 0(/3 2 ): 

/3 2 8 2 
E ^ eo + ei tL F ^ /o + /l ^_ (5.5) 

with 

eo = cto + a i an d /o = a o ~~ °>i- (5-6) 

Since K q = K c + 5K C , we deduce that 

^ = -^sinh^). (5.7) 

Thence, expanding to 0(e), we find, 

1 1 ei ei 



cosh # + 1 cosh 9 cosh 9 + sin ai7r cosh 9 — sin ai7r 
eiecosai7r ( 1 1 



T{9) 



2sinax7r \ (cosh 9 + sin aiii) 2 (cosh 9 — sin ai7r) 2 

+ . 6/1 (5-8) 
smai7rcosn 9 ) 

Comparing ( |5.2|) with (|5.8| ) we see a pleasing similarity. In fact the two formulae 
are identical, to 0(e), provided we choose e\ and j\ suitably. In other words, we may 
deduce that 

e 1 = -2a 1 + — = -(a + ai) +0(e 2 ), (5.9) 

7rsmai7r 

and fx is proportional to e. Unfortunately, the calculation does not allow anything 
more detailed to be learned about f\. To do better would need a correction to the 
reflection factor to 0(e 2 ). 
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6 Discussion 



The purpose of this calculation was to test a little more deeply the expression for the 
sinh-Gordon particle reflection factor given in || and to learn additional information 
concerning its dependence on the boundary parameters a and <j\. The result of the 
investigation is gratifying because it agrees with alternative derivations of the reflection 
factor and it also agrees with the following conjecture. Everything we have learned so 
far is consistent with quite simple expressions for E and F: 

E = (a + a 1 )(l-B/2) F = (a - a x )(l - B/2), (6.1) 



where the coupling constant dependence enters via the expression for B given in 
Similar expressions for these parameters have been arrived at via other arguments by 
Zamolodchikov |3"3fl . 

If ( |6.1| ) is correct then the reflection factor is invariant under the interchange ao <-> 
a\. In effect, this invariance restores the Z2 bulk symmetry which apparently was 
broken by the boundary condition and replaced by a symmetry under the simultaneous 
interchange of </> with —cf) and a with a\. The reflection factor is also invariant if ao 
and/or a>i is replaced by its negative, as it should be given the definitions of a and 
<j\. It is consistent with what is known at the special value of the coupling constant, 
known as the 'free-fermion' point in the sine-Gordon model, where B = —2 and the 
S-matrix is unity. There, the restrictions on the parameters in the reflection factor can 
be solved exactly and are in agreement with (|6.1| ) f23 . 



Note also, that with the expressions ( |6.1| ) the reflection factor Q3.2| ) has a weak- 
strong coupling symmetry matching the symmetry of the S-matrix under (3 — > 47r//3. 
To see this, note that setting 

4:71 

(a* ,a* 1 ,P*) = ^(a ,a 1 ,(3) (6.2) 

defines a new triple of coupling constants with the property that 

K q {6, a , a 1; (3) = K q {6, a*, a*, f3*). (6.3) 



If ( |6.1| ) is correct, implying the duality symmetry ( |6.2| ), then we are faced with other 
puzzles. For example, it is known that the supersymmetric version of the sinh-Gordon 
model is only integrable when restricted to a half-line with some very special boundary 



conditions (either a = ai = or, a = ai = 1) (see f34j), and this restriction would 



appear to be incompatible with a weak-strong coupling symmetry without modifying 



It is also known HQ that the other affine Toda field theories constructed from data 
in the ade series, when restricted to a half-line, allow only a finite number of possible 
boundary conditions. In fact, the a^ or sinh-Gordon model is apparently the only 
example within this series which allows continuous boundary parameters. Expressions 
for the associated reflection factors for the other models are largely unknown but, it 
will be interesting to discover if they too can permit a duality symmetry in the presence 
of a boundary which will match the symmetry of their bulk S-matrices. 
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